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Abstract 

Using the Painleve analysis preceded by appropriate transforma- 
tions of nonlinear systems under investigation, we discover two new 
cases in which the Pietrzyk-Kanattsikov-Bandelow vector short pulse 
equation must be integrable due to the results of the Painleve test. 
Those cases are technologically important because they correspond to 
the propagation of polarized ultra-short light pulses in usual isotropic 
silica optical fibers. 



The short pulse equation (SPE), which has the form 

u,t = u + l{u')^^ (1) 

up to a scale transformation of its variables, was introduced recently by 
Schafer and Wayne [1] as a model equation describing the propagation of 
ultra-short light pulses in silica optical fibers (note, however, that for the first 
time this equation appeared in differential geometry, as one of Rabelo's equa- 
tions describing pseudospherical surfaces p]). Unlike the celebrated nonlin- 
ear Schrodinger equation which models the evolution of slowly varying wave 
trains, the SPE is well applicable when the pulse spectrum is not narrowly 
localized around the carrier frequency, that is when the pulse is as short as 
a few cycles of its central frequency. Such ultra-short pulses are very impor- 
tant for future technologies of ultra-fast optical transmission of information. 
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The SPE dl]) is an integrable equation possessing a Lax pair [21 S] of the 
Wadati-Konno-Ichikawa type [5]. The transformation between the SPE and 
the sine-Gordon equation was discovered in |3], and later it was used in [6] 
for obtaining exact loop and pulse solutions of the SPE from the well-known 
kink and breather solutions of the sine-Gordon equation. The derivation 
of that transformation was considerably simplified in [7], where analogous 
transformations to well-studied equations were also found for the other three 
equations of Rabelo. The recursion operator [1], Hamiltonian structures and 
conserved quantities [HI [9], Hirota's bihnear representation [IQl, multisoliton 
solutions [11] and periodic solutions [12], [13] of the SPE ([1]) were found and 
studied as well. 

Very recently, Pietrzyk, Kanattsikov and Bandelow [2] introduced the 
vector short pulse equation (VSPE), a two-component nonlinear wave equa- 
tion that generalizes the scalar SPE ([1]) and describes the propagation of 
polarized ultra-short light pulses in cubically nonlinear anisotropic optical 
fibers, which can be written as 

Un,xt = CniUi + Cnijk (UiUjUk)^^ , (2) 

where n,i,j,k = 1,2, and the summation over the repeated indices is as- 
sumed. Since the constant coefficients c„i and Cnijk are determined by optical 
properties of the fiber's material, there is a wide variety of mathematically 
different cases of the VSPE and it is interesting to find out which of them 
are integrable systems of coupled nonlinear wave equations. Three integrable 
cases of the VSPE (|2]) were obtained in [H] by direct construction of their 
Lax pairs, namely, 

u^t = u + I {u^ + 3uv^) , v^t = V + I {3u'^v + v^) , (3) 
u^t = u + I {u^ - 3uv^) , v^t = V + I {3u'^v - v^) , (4) 

and 

where u and v denote the polarization components Ui and U2- The following 
valuable remark was made in [14J on the nature of the system ([5]): this case 
of the VSPE describes the propagation of a small perturbation v on the 
background of a solution u of the scalar SPE ([1]). Contrary to what was 
proposed in |14.|, however, we cannot consider the VSPE ([3]) as a short pulse 
counterpart of the Manakov system [15] of coupled nonlinear Schrodinger 
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equations. Indeed, in the new variables p = u + v and q = u — v the equations 
of the system ([3D become uncoupled and turn into two scalar SPEs ([T]) for 
p and q separately, whereas the polarization modes in the Manakov system 
do interact nonlinearly [16] . In the variables p = u + iv and q = u — iv the 
equations of the system (jl]) become uncoupled as well. It is easy to prove 
that the 4x4 zero-curvature representations, found in [T3] for ([3]) and (jl]), 
can be brought by gauge transformations into the block-diagonal form with 
the 2x2 diagonal blocks corresponding to the zero-curvature representation 
of the scalar SPE ([1]) for p or q. 

In the present paper, we show that there are at least two more cases of 
the VSPE 02]) which can be strongly expected to be integrable due to the 
analytic properties of their general solutions, namely, 

u^t = u + \[u^ + uv'^)^^, v^t = v + \[u^v + v^)^^, (6) 

and 

Uxt = u + \ {u^) , v,t = v + \ . (7) 

The VSPE ([6]) represents the technologically important case, where the fiber 
is made of a cubically nonlinear isotropic optical material, such as the widely 
used silica glass, but the ultra-short light pulse is not linearly polarized; if 
integrable, this system can be interesting as a short pulse counterpart of the 
Manakov system. The VSPE ([71) can be considered as the limiting case of the 
system ([6]) for small values of f , that is the case of almost linearly polarized 
pulses. We discover these systems ([6]) and ([7j) by applying the Painleve test 
for integrability of partial differential equations [TTl [HI [19] to the following 
two one-parameter classes of VSPEs: 

u^t = u + \ [u^ + cuv"^) , v^t = v + \{cu\ + v^)^^, (8) 

and 

«xt = M + I v,t = v + l[u'v)^^, (9) 

where c is the parameter. Of course, this is far not the complete test for 
integrability of the whole variety of systems ([2|) but rather the first attempt 
to search for new integrable VSPEs systematically. In what follows, we make 
the computations using the Mathematica computer algebra system |20] and 
omit their inessential bulky details. 

Let us consider the class of systems ([8]) first. It is easy to see that the 
Painleve test cannot be applied to the VSPE ([S]) directly, for the reason of 
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an inappropriate dominant behavior of solutions near a movable singularity 
manifold, and we must appropriately transform the nonlinear system under 
investigation in order to improve this behavior and start the test. This is a 
point of crucial importance in our study. We follow the way of transformation 
similar to the way used in [7j for the scalar SPE ([1]) and other Rabelo's 
equations. We make the change of the independent variable x, 

uix, t) = f{y, t), v{x, t) = g{y, t), y = y{x, t), (10) 

and determine the function y{x,t) by the relation 

yt = l{u^ + v^)y,. (11) 

In this relation flTTl) . the polynomial u"^ + f ^ is taken for the reason of sym- 
metry between u and v; one could use there a general quadratic polynomial 
in u and v instead, but this would have no effect on the dominant behavior 
of solutions, positions of resonances and compatibility of recursion relations, 
found during the Painleve analysis. Then, inverting y = y{x,t) as x = x{y, t), 
we obtain from (jS]) and (ITT!) the following system of three coupled equations 
for fiy,t), g{y,t) and x{y,t): 

2xt + f + g^ = 0, 

Q^lfyt + (3 - C)g^Xyfyy " 2cfgXygyy + ((C " 3) g^ fy + 2cfggy) Xyy 

+ (6 - 4:c)gfygyXy - 2cfglxy - = 0, 

^Xlgyt + (3 - C)fXygyy " 2CgfXyfyy + {{C~ 3) f gy + 2Cg f fy) Xyy 

+ (6 - 4c)fgyfyXy " 20^/2^^ - 6gxl = 0. (12) 

Note that the fact of correspondence between the fifth-order system f lT^ and 
the fourth-order system (we mean here the total order of a system, or the 
number of arbitrary functions in its general solution) can be explained by 
the invariance of (fT2|) with respect to an arbitrary transformation y \^Y{y), 
which just means that the solutions of (fT2|) represent the solutions of ([8]) 
parametrically, with y being the parameter. 
Substitution of the expansions 

X = xo{y, t)<p{y, tr + --- + xr{y, t)<j){y, ty+'' + ■ ■ • , 
/ = /o(i/, tf + --- + friy, my, ty^P + ■ ■ ■ , 
g = go{y, mv, ty + ■ ■ ■ + g^y, mv, ^Y^^ + ■■■ (13) 
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to the system fll2p determines the dominant behavior of solutions in the 
neighborhood of a manifold (p{y,t) = 0, i.e. the admissible values of a, jS, 
7, xo, fo and go, and the corresponding positions of resonances r, where 
some arbitrary functions can enter the expansions. If we assume that at 
least one of the exponents a, /3, 7 is negative, we immediately get a = 
/3 = 7 = — 1 for all values of c except c = — 1 (we do not know at present 
how to transform the systems ([H]) and (jH]) with c = — 1 in order to start 
the Painleve test for them, and we do not consider the case of c = — 1 in 
this paper). We have to exclude the manifolds = with 0^0^ = 0, for 
which no well-posed recursion relations appear for the coefficients x„, /„, gn 
of the expansions (IT^ : the reason to exclude such characteristic manifolds 
consists in that arbitrarily nasty singularities of solutions can occur along 
characteristics, even for integrable equations [211 122]- For non-characteristic 
manifolds, without loss of generality, we choose (j)y{y,t) = 1 with 0t 7^ and 
set all the coefficients Xn, fn, gn in ( fT3|) to be functions of t only. Then we 
find for c 7^ —1, 1, 3 that 

Xo = -(1 + c)0t, fo = ±Wl + c0t, go = ±iVl + c4>t, 

r = -1, 1, 4, i (5 - yi^) , I (5 + yi^) , (14) 

where the ± signs in the expressions for fo and go are independent, and 
i^ = —1. The case of c = 3, i.e. the VSPE ([3]), is not of interest because the 
equations can be easily uncoupled. And for c = 1 we find that 

xo = -2(j)t, fl + gl = -AcPl r = -1,0,1,4,5, (15) 

where goif) (or foif)) is arbitrary due to the resonance r = 0. In the case 
of relations ([SD, there are only three possibilities to have four resonances 
in integer non-negative positions, namely, c = —1,0, 1. However, the values 
c = —1,1 have been excluded, whereas the value c = corresponds to the 
uninteresting case of uncoupled equations in the VSPE (|8]). On the other 
hand, in the case of relations fll5p corresponding to c = 1, we have already 
got the admissible positions of resonances. Continuing to study this case, 
we derive from the equations f|T2|) with c = 1 the recursion relations for the 
coefficients of expansions f|T3l) with a = /3 = 7 = — 1, check the compatibility 
of those relations at the resonances, and find that all the compatibility con- 
ditions are satisfied identically. Thus, in this case, the singular expansions 
of solutions turn out to be some generalized Laurent series containing no 
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non-integer powers and no logarithmic terms, the arbitrary functions in the 
series being go{t), Xi{t), g^{t), g^{t), and ^/'(t) in = |/ + ^(t). 

To complete the Painleve test for the system f|T2|) with c = 1, we also 
have to study the expansions f|T3|) with a = /5 = 7 = 0. Such expansions, 
which start like Taylor series, exist for every partial differential equation. 
Usually, they are the Taylor expansions of regular solutions governed by the 
Cauchy-Kovalevskaya theorem [23] • In some cases, however, such Taylor- 
like expansions can contain non-dominant singularities; this may happen, 
for example, when the Kovalevskaya form of a studied equation is singular 
for some of Cauchy data [211 l25|. Substituting the expansions (fT3|) with 
a = /? = 7 = 0to the system (IT^ with c = 1, we find that no well-posed 
recursion relations appear for the coefficients x„, /„, Qn if the manifold = 
is determined by the conditions (pyCpt = or (/q +5^0 + 6a;o,t) 0?/ = 6xo^y(j)f 
We exclude such characteristic manifolds, choosing (py = 1 with 0* 7^ and 
setting all the coefficients a;„, /„, gn in f|T3|) to be functions of t only, with 
/o + S'o + 63^0,* 7^ 0. Then we find that the positions of resonances depend 
on whether the Cauchy data satisfy the condition /q + (^q + '^^o,t = 0. When 
/o + S'o 7^ "2x0,4, the arbitrary functions in the expansions are XQ{t), fo{t), 
doit), fi{t), gi(t), as well as ■0(t) in = y + ip(t) (the appearance of one 
extra arbitrary function in Taylor-like expansions of solutions was discussed 
in [23]). When /q + (^o ~ ~2xo,f, the arbitrary functions in the expansions 
are Xo(t), go(t), gi(t), g2(t), as well as tp(t) in (p = y + ip{t). In both cases, 
the positions of resonances are integer, and the recursion relations are com- 
patible at the resonances. Therefore the expansions are some Taylor series 
containing no non-integer powers and no logarithmic terms. Consequently, 
the system f|T2l) with c = 1 has passed the Painleve test, and the VSPE ([6]) 
can be strongly expected to be integrable. 

Let us proceed now to the class of systems Making the same trans- 
formation (fTOl) . setting 

yt = (16) 

and inverting y = y{x,t) as x = x{y,t), we obtain from ([9]) and f|T6l) the 
following system of three coupled equations for f{y,t), g{y,t) and x{y,t): 

2xt + f = 0, fyt - fXy = 0, 

^^igyt + (3 - C)fXygyy " 2CgfXyfyy + - 3)fgy + 2Cgffy) Xyy 

+ (6 - ^c)fgyfyXy " 2cgf^Xy - 6gxl = 0. (17) 
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Then, using the expansions (fT5]) . choosing (py = 1 with 0^ 7^ to exclude 
characteristic manifolds = 0, setting all the coefficients x„, /„, Qn to be 
functions of t only, and assuming that at least one of the exponents a, jS, 7 
is negative, we find from the system f|T71) that 

a = (3 = -l, 7^-87 + 2 = 6/0, 
xo = -20t, /o = ±2i0t, r = -1,0,1,4,3-27, (18) 

where the resonance r = corresponds to the arbitrariness of go{t). The 
resonance r = 3 — 27 must correspond to the arbitrary coefficient g^^2'y{t), 
due to the structure of recursion relations which follow from (fT7|) and (fT3|) . 
Denoting the resonance position 3 — 27 as m and using (fT8|) . we have 

7 = (3-m)/2, c = 24/(m2-l). (19) 

The admissible values of m are m = 2, 3, 4, 5, ... , since m = is excluded 
because it leads to the double resonance r = 0, which indicates that the ex- 
pansion for g must contain a logarithmic term, and m = 1 is excluded because 
it implies c = 00. Though the even values of m correspond to non- integer 
values of 7, these cases should not be excluded, because by introducing the 
new variable h = g'^ one can improve the dominant behavior of solutions. 

We have found infinitely many cases of the system f[T7|) . which are all char- 
acterized by some admissible dominant behavior of solutions and admissible 
positions of resonances. Unfortunately, we cannot check the compatibility of 
the recursion relations at the resonances for the whole infinite set of those 
cases at once. This situation is quite similar to the one observed in the 
Painleve analysis of triangular systems of coupled Korteweg-de Vries equa- 
tions On available computers, we were able to complete the Painleve 
test for the cases m = 2, 3, . . . , 9, 10 of the system f|T7j) with c given by f|T9|) . 
The recursion relations turn out to be compatible only in the cases m = 3 
and m = 5, whereas some nontrivial compatibility conditions appear in all 
other cases at the resonance r = m as an indication of non-dominant log- 
arithmic singularities of solutions. The case m = 3 with 7 = and c = 3 
corresponds to the integrable VSPE discovered in [H]. The case m = 5 
with 7 = — 1 and c = 1 corresponds to our new VSPE ([7]). Continuing 
to study the system (fT7|) with c = 1, we consider the expansions (fT3|) with 
a = /? = 7 = 0, and they turn out to be some Taylor series containing no 
non-integer powers and no logarithmic terms. Consequently, the system (1171) 
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with c = 1 has passed the Painleve test, and the VSPE ([7j) can be strongly 
expected to be integrable. 

Let us remind, however, that the Painleve property does not prove the 
integrability of a nonlinear equation but only gives a strong indication that 
the equation must be integrable. Consequently, the new probably integrable 
nonlinear systems and ([7]), discovered in this paper, deserve further inves- 
tigation, especially taking into account their importance for physics and tech- 
nology. To find their Lax representations, generalized symmetries, Hamilto- 
nian structures and soliton solutions seems to be a complicated problem, 
but a sufficiently interesting one to attract attention of experts in nonlinear 
mathematical physics. 

The main part of this research was carried out during the author's visit 
to the Max Planck Institute for Mathematics (Bonn, Germany), whose hos- 
pitality and support are acknowledged with a deep gratitude. 
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